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Abstract. Mass and radius of planets transiting their host stars are provided by radial ve- 
locity and photometric observations. Structural models of solid exoplanet interiors are then 
constructed by using equations of state for the radial density distribution, which are compliant 
, with the thermodynamics of the high-pressure limit. However, to some extent those structural 
models suffer from inherent degeneracy or non-uniqueness problems owing to a principal lack of 
knowledge of the internal differentiation state and/or the possible presence of an optically thick 
^ , atmosphere. We here discuss the role of corresponding measurement errors, which adversely 
affect determinations of a planet's mean density and bulk chemical composition. Precise mea- 
surements of planet radii will become increasingly important as key observational constraints 
for radial density models of individual solid low-mass exoplanets or super-Earths. 
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. 1. Introduction 

The steadily growing number of transiting planet discoveries will allow us to charac- 
terize rocky exoplanets in terms of internal structure and atmospheric composition with 
important implications for their formation, orbital evolution, and possible habitability 
(see review by Haghighipour 2011). Numerical models of planetary interiors based on 
laboratory data of physical material properties are aimed at improving the general un- 
derstanding of their origins, pathways of evolution, and diversity. In case of the terrestrial 
planets and satellites within the solar system, the resultant radial profiles of density and 
related material properties are required to satisfy geophysical observations and cosmo- 
chemical arguments for the compositions of major geochemical reservoirs such as the 
core, mantle, and crust (see Sohl & Schubert 2007, Sohl 2010, and references therein). 
^{■^ , For rocky exoplanets, the numerical models have to be consistent with the observed plan- 
ar - ) 1 etary masses and radii measured from ground-based observations and space missions. 

Calculated models have been used to derive mass-radius relationships for exoplanets as- 
suming a range of chemical compositions to gain insight into the bulk compositions and 
possible interior structures of these planets (e.g., Valencia 2011, and references therein). 



Wagner et al. (2011a) have reinvestigated mass-radius relationships for rocky exoplanets 
using equations of state that are compliant with the thermodynamics of the high-pressure 
limit of a given material. 

The mean density is the main indicator of the bulk composition of solid planets. The 
purpose of this paper is to discuss the role of mass and radius measurement errors for 
determinations of a planet's mean density and bulk chemical composition by using cal- 
culated relationships between radius and mass of solid low-mass exoplanets. 



2. Method 

Since there are usually fewer constraints than unknowns, even basic interior structure 
models that would involve only two or three chemically homogeneous layers of constant 
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density suffer from inherent degeneracy. In addition, structural models of low-mass ex- 
oplanets have non-uniqueness problems because of their unknown differentiation state, 
and/or the possible presence of an optically thick atmosphere (Valencia 2011), and the 
extrapolation of equations of state of mineral phase assemblages to extremely high pres- 
sures (Swift et al. 2012). Further uncertainties are related to pressure-induced mantle 
phase transitions, the stability field of post-perovskite, the possible presence of addi- 
tional high-pressure silicate phases, and the physical state of core iron alloys (Fortney 
et al. 2009). This has important implications for mass-radius relationships and their us- 
age for the characterization of low-mass exoplanets in terms of bulk composition and 
interior structure and the possible existence of self-sustained magnetic fields (Wagner 
et al. 2011b). 

To obtain relationships between planetary radius R p and total mass M p for solid 
planets according to 

R P = cM?, (2.1) 

where j3 is the scaling exponent, we hrst solve the structural equations for mass and 
momentum assuming hydrostatic equilibrium (e.g., Wagner et al. 2012) 

dm o . . dP Gm , . . . 

— - 47rr 2 p(r); — = r p(r), (2.2) 

dr dr r z 

where G is the gravitational constant, together with an equation of state (EoS) for the 
local density p according to 

p(r) = f EoS (P). (2.3) 

Here, we adopt the generalized Rydberg EoS with material parameters that are valid 
for iron (Fe), silicate (MgSiOa), and water ice (H2O). This EoS has the advantage of 
being compliant with the high-pressure limit of a given material, thereby describing 
the radial distribution of density p(r) up to extremely large compression ratios of solid 
planetary materials (see Wagner et al. 2011a, and references therein). We further neglect 
the thermal pressure contribution because of its minor importance for massive exoplanet 
interiors (e.g., Seager et al. 2007). Successful solutions of the structural equations are 
required to satisfy the boundary conditions 

m(0)=0, P(Rp)=0, P(r pt )=f P t (2.4) 

at the planet's center (r = 0), surface (r = R p ), and pressure-induced phase transition 
boundaries (r = r p t). We are taking into account transition pressures f p t for the a- and 
e-phase of iron, the perovskite and post-perovskite phase of MgSiC>3, and the low- and 
high-pressure phases of water ice. 

Table 1. Scaling exponent ft for solid exoplanets in the mass range of 1 to 10 M®. 



Composition 


|Valencia et al. 


(2006)| |Sotin et al. (2007)| |Wagner et al. (2011a)| 


Earth-like 


0.262 


0.274 0.267 


Ocean planet 


0.244 


0.275 0.261 


Mercury-type 


none 


none 0.269 



Note: Assumed is a radius-mass relationship according to R p — cM^ with scaling exponent /?. 

Second, a power law fit is performed to adjust the scaling exponent /3 in eq. 12.11 for 
solid low-mass exoplanets in the mass range of 1 to 10 Mg. Tab. [T] provides a comparison 
of published scaling law exponents for solid low-mass exoplanets of various compositions. 
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Figure 1. Mass-radius diagram for planets with different bulk compositions compared to cur- 
rently known low-mass exoplanets in Earth units. We divide equilibrium surface temperatures 
into three domains from 500 to 1000 K; 1000 to 1500 K; and 1500 to 2000 K. While the solid 
curves denote homogeneous, self-compressible solid spheres of water ice, silicate rock, and iron, 
respectively, the dashed curves exhibit differentiated models of intermediate bulk compositions. 

Because of the large compression ratios involved, those are substantially less than (3 = 1/3 
in case of a homogeneous density distribution. 

The mean density of a spherical planet is given by 



3 M p 



(2.5) 



where, in general, mass M p and radius R p of planets transiting their host stars are 
provided independently from each other by radial velocity and photometric observations. 
We therefore employ an error propagation analysis according to 



Ap = 



dp 



AM P 2 



dp_ 

8R rj 



(2.6) 



and obtain, upon substitution of the radius-mass relationship given in eq. 12.11 the prop- 
agated relative error in mean density 



(2.7) 



as a function of the scaling law exponent (3 and the observational uncertainties of the mass 
and radius determinations. The latter can be expressed in terms of some key observables, 
namely 
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where K* and S* denote radial velocity semi-amplitude and transit depth of the host star, 
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Figure 2. Relative error in mean density, Ap/p, as a function of relative uncertainties in plan- 
etary mass, AMp/M p , and radius, AR P /R P , as given by eq. 12.71 Assumed is a radius-mass 
relationship according to R p = cM^ with scaling exponent /3 = 0.27, representative for solid 
low-mass exoplanets or super-Earths. 

respectively. Upon insertion of (3 — 1/3 in eq. 12. 71 for a homogeneous density distribution, 
the relative error in mean density can be written as 



Ap 
P 



/ AMp 
V M p 



9 



(2.9) 



which is a valid approximation for small terrestrial planets and moons (M p 1 M ffl ) 
that are subject to low internal pressures and negligible compression ratios. 



3. Results and Discussion 

The mass-radius diagram shown in Fig. [T] exhibits the remarkable compositional di- 
versity of known transiting low-mass exoplanets. Whereas the low density of GJ 1214b, 
that has no analogue in the solar system, is attributed to the presence of an extended 
optically thick hydrogen atmosphere (Rogers & Seager 2011), it has been suggested that 
carbon-rich solids may predominate the bulk composition of 55 Cancri e (Madhusudhan 
et at 2012). CoRoT-7b and Kepler-lOb, on the other hand, are thought to be composed 
of silicate rock and iron with the latter concentrated in massive metallic cores (Wagner 
et at 2012, Swift et al. 2012). These compositional trends are roughly followed by the 
calculated equilibrium surface temperatures, with relatively hot planetary environments 
along and below the silicate line and colder but still hot environments along and above 
the water ice line. Furthermore, it is important to note that solid planets with masses 
above 10 M ffi were not discovered to date, although transit detections should be readily 
feasible owing to the relatively large planet radius. 



Mass-radius relationships of rocky exoplanets 



■5 



35,000 




generalized Rydberg 
Keart9 
reciprocal K' 



0.5 1 

Radius, R/R 



0.5 1 

Radius, R/FV. 



Figure 3. Radial distributions of (a) density and (b) gravitational acceleration for solid model 
planets of 5 Mq (lower set of model curves) and 10 M® (upper set of model curves), respectively. 
Planet bulk composition is held fixed to that of the Earth. Using the generalized Rydberg, Keane, 
and reciprocal K' EoS, the resultant surface radii differ by less than 2% from each other (adopted 
from Wagner et al. 2011a). 



Inspection of Tab. Q] suggests that the scaling exponent f3 is similar for radius- mass re- 
lationships of solid exoplanets in the mass range of 1 to 10 M®, irrespective of their bulk 
compositions. To calculate ellipsoidal distributions of the relative error in mean density, 
Ap/p, we thereby assume an average scaling exponent ft = 0.27, being representative for 
solid low-mass exoplanets or super-Earths. From inspection of Fig. [21 it is seen that the 
accurate determination of mean planetary density depends on the observational uncer- 
tainties of both planetary mass and radius. The performance comparison of ongoing and 
proposed transit missions given in Tab. [2] hints at the expected accuracy of correspond- 
ing planet density determinations. Owing to the variable slope of mass-radius relations, 
planetary mean density in the giant-planet mass range is mainly constrained by precise 
radius determinations, whereas mass and radius pose equally important constraints on 
structural models of low-mass planets. 



Table 2. Performance of ongoing and proposed transit missions. 

Uncertainty CoRoT a Kepler 6 CHEOPS c PLATO d 

radius, ±0.06 ±0.02 ±(0.02-0.1) ±0.02 

mass, ±0.15 ±0.27 ±0.10 ±0.10 

" (CoRoT-7b) (Kepler-lOb) 



References: 

a) |Bruntt et a l. (2010)] b) |Batalha et al. (2011)1 c) CHEOPS proposal Qhttp : //cheops . unibeTch) ; d) PLATO 
M3 proposal (http : //www . oact . inaf . it/plato/PPLC/Home .html) . 



In Fig. [3j structural model determinations of atmosphere-less solid planets are com- 
pared for different equations of state and fixed Earth-like bulk compositions. It is seen 
that the trade-off in calculated planet radius is smaller than present-day measurement 
uncertainties from transit photometry, indicating that mass-radius relationships of solid 
exoplanets are sufficiently robust to infer their bulk compositions (Wagner et al. 2011a). 
Since mass determinations using the radial velocity method are currently limited to 
an accuracy of about ±10 percent, precise measurements of planet radius will become 
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increasingly important as key observational constraint for radial density models of in- 
dividual low-mass exoplanets. At any rate, the precise and reliable characterization of 
planet host stars in terms of stellar radius, mass, and age is prerequisite to infer the 
mean density and thereby composition of low-mass exoplanets or super-Earths with con- 
fidence. This can be accomplished by using asteroseismology on sufficiently bright and 
small target stars (e.g., Rauer et al. 2011). 

4. Conclusions 

Current detection limits of ground-based observational methods have limited the dis- 
covery of solid exoplanets to only a few, although low-mass planets beyond the solar 
system should be quite abundant (e.g., Marcy et al. 2011). Mass-radius relationships 
based on numerical models of solid exoplanet interiors are sufficiently robust to infer 
a planet's bulk composition from accurate determinations of its mean density. The fu- 
ture detection of transiting mini-Neptunes and super-Earths will provide fundamental 
information to better constrain their bulk compositions and possible interior structures 
in terms of metallic cores, silicate/carbon-rich mantles, water-ice/liquid shells, and/or 
atmosphere mass fractions. 
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